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Abstract We describe a recently realized experiment 
producing the most spherical cavitation bubbles today. 
The bubbles grow inside a liquid from a point-plasma 
generated by a nanosecond laser pulse. Unlike in previ- 
ous studies, the laser is focussed by a parabolic mirror, 
resulting in a plasma of unprecedented symmetry. The 
ensuing bubbles are sufficiently spherical that the hy- 
drostatic pressure gradient caused by gravity becomes 
the dominant source of asymmetry in the collapse and 
rebound of the cavitation bubbles. To avoid this natu- 
ral source of asymmetry, the whole experiment is there- 
fore performed in microgravity conditions (ESA, 53rd 
and 56th parabolic flight campaign). Cavitation bub- 
bles were observed in microgravity (~0g), where their 
collapse and rebound remain spherical, and in normal 
gravity (Ig) to hyper-gravity (1.8(7), where a gravity- 
driven jet appears. Here, we describe the experimental 
setup and technical results, and overview the science 
data. A selection of high-quality shadowgraphy movies 
and time-resolved pressure data is published online. 
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1 Introduction: A roadmap towards unification 

Cavitation bubbles remain a key topic in fluid dynam- 
ics. While traditionally associated with erosion dam- 
age ( Philipp and Lauterbornl . 1998]), they are currently 
reconsidered in a wide range of modern applications 
wit hin food technology ( Mason et all 19961 ). water clean- 
ing ( Wolfrum et all. 2003 ). medicine (Leighton and Clev e 



land, l2010|), and microfluidics (jTandiono et al 1201 lh . 



This extraordinary breadth of applications originates 
from the rich physics governing the collapse of indi- 
vidual bubbles. Fundamental science is now asked to 
provide a robust framework unifying this richness. 

All research on cavitation has departed from the 
ideal model of a perfectly spherical bubble collapsing 
within a liquid medium. The first time-solutions for 
the case of an empty bubble in an infinite , invis cid, 
incompressible liquid dStokesl Il847l : iRavleighl . Il917l ) ex- 
hibit infinite velocities of the bubble wall at the collapse 
point . This branch point singularity (jObreschkow et al 
20121 ) highlights the insufficiency of the Stokes- Rayleigh 
approach, and resulted in ever more detailed consider- 
ations of the complex phenomena triggered during the 
bubble collapse. Those post-collapse processes can be 
grouped into four classes (see Fig. [TJ: 

(a) rebound bubbles a rising when bubbles bounce off 
their enclosed gas (jAkhatov et al 1200 1[ ) ; 

(b) micro-jets emerging in an asymmetric collapse (e.g. , 
Blake et al 1 19991 : Katzl . 1 19991: I Wang and Blakel . 



20ld : IObreschkow et all . 120061 : feobel et alil2009f h 
(c) shock waves caused by the liquid c ompr ession at 
the collapse point (e.g.. Ohl et all. 1 19991: Schnerr 
et al. 120081 : lObreschkow et all . l2011al ): 
thermal effects, e.g. , heating/cool i ng (Hickling| . ll965l ). 
luminescence ( e.g. . iBrenner et all . I2002T ) 
ical reactions (|Suslickl . Il990l ). 



(d) 



and chem- 
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Fig. 1 Illustration of the four classes of phenomena triggered by collapsing cavitation bubbles. All images were obtained using 
the setup presented in this paper, but they correspond to different experimental conditions. Not all phenomena are active in 
every collapse. For example, a spherical collapse in microgravity produces no jets. 



In the wealth of established facts on those post- 
collap se phenomena (see review by Lauterborn and Kurd . 
2010), we might still be missing a unifying picture. Most 



prominently, a general theory for the relative impor- 
tance of the different phenomena in various experimen- 
tal conditions remains to be uncovered. How does the 
energy partition between rebound, shock, jet, and lu- 
minescence vary as a function of the liquid pressure, 
the liquid temperature, the gas content and the bub- 
ble sphericity? Quantitative answers to these questions 
promise to become an outstanding tool for optimizing 
virtually any application that relies on a particular fea- 
ture of collapsing cavitation bubbles. 

On the experimental level, three requirements must 
be met in the quest for a unified theory of the energy 
relocation of collapsing bubbles: 



This paper presents an experiment meeting those 
three criteria, and it provides online access to a se- 
lection of experimental data. Section [2] describes the 
experimental setup, its peculiar features, and the con- 
ditions in which the experiment was performed. Sec- 
tion |31 then presents results that quantify the quality 
of the data and illustrate the wealth of observed phe- 
nomena. For example, we present high-speed visual- 
izations of the gravity-driven jet produced during the 
collapse of a spherical cavitation bubble in water sub- 
jected to normal gravity. Systematic scientific results 
regarding the jet formation and the energy partition 
between rebo und and shock were presented in separate 
manu scripts ([Obreschkow et all . l2011bl ; iTinguelv et all . 



20121 ). Section @] explains the data structure and access 



to the online data. Section [S] concludes the paper. 



The setup must allow the generation of highly spher- 
ical bubbles that conserve their sphericity in the re- 
bound, hence suppressing any jets. Jet formation 
can then be stimulated by adding controlled asym- 
metries. Only through such a precise control can 
different energy channels be disentangled. 
The setup must be equipped with sensors measur- 
ing the essential phenomena: a high-speed camera 
to capture rebound bubbles and potential jets, a 
color-sensitive light sensor to measure the intensity 
and temperature of sonoluminescent radiation, and 
a time-resolved pressure sensor to capture shocks. 
The experiment must be performed for a wide array 
of initial conditions to sample the parameter space 
covered by typical applications. 



2 Experimental setup 

The experimental setup can generate a spherical va- 
por bubble (radius = — 7 mm) in the middle of a cu- 
bic liquid volume (178 x 178 x 150 mm 3 ). The growth, 
collapse, luminescence, and rebound of this bubble are 
filmed using a high-speed camera, while shock waves 
are recorded by a pressure sensor. To modulate and 
annihilate the gravity-induced pressure gradient in the 
liquid, the experiment is performed inside an aircraft 
performing parabolic flights. Section POl overviews the 
mechanical structure of the setup, whose subsystems 
are described in Sections l2.2H2.5l The flight manoeuvres 
are explained in Section [521 and the time-sequence of a 
single experiment is detailed in Section |2~71 Additional 
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background i nformation on the experimental setup is 
provided by ( Tinguelv . 20131 ). 



2.1 Structure of the experimental setup 

Fig. [2] illustrates the experimental setup. The entire 
setup is contained within a mechanical structure called 
"the rack" (see Fig. [5^). Its robust skeleton of strut pro- 
files (50 x 50 mm 2 and 30 x 30 mm 2 ) and two horizon- 
tal aluminium plates (thickness 10 mm) ensures that 
precision-parts, such as optical components, displace 
less than 10 /im, when passing from normal gravity 
to weightlessness and hypergravity, respectively. Ulti- 
mately, the whole rack withstands accelerations up to 
9g as required for flight security. The rack is fixed inside 
the aircraft "A300 zero-g" close to the aircraft's center 
of gravity (Fig. where the best level of micrograv- 
ity can be achieved (see Section |2~B|) . Globally, the rack 
consists of two parts: an open "lower stage" and a closed 
"upper stage" , which can be opened via a lid. A laptop 
is fixed on top of this lid. A custom-designed Lab View 
script running on this computer controls the experimen- 
tal conditions, records the sensor data, and automati- 
cally triggers the bubble generation at the preselected 
level of gravity. 

The upper stage of the rack is schematically dis- 
played in Fig. [5J: and detailed in Figs. [5b, d. This stage 
contains three main systems: a pressure-controlled test 
chamber filled with liquid (Section 12. 2j) , a laser-system 
to generate a cavitation bubble inside the test cham- 
ber (Section [23]), and a high-speed imaging system to 
record the evolution of the cavitation bubble (Section 
12. 4[) . Additionally, the upper stage also hosts various 
additional sensors, described in Section [2.51 

The lower stage of the rack hosts auxiliary electron- 
ics and power supplies. A USB-interface (Phidgetlnter- 
faceKit 8/8/8 #1018) converts analog sensor signals 
(Section 12. 5[) into computer-readable digital data. This 
interface also permits the computer to release a laser 
pulse for bubble generation (Section l2.3[) . to control the 
pressure pump (Section 12. 2[) . and to activate the illu- 
mination system (Section 12.41) . An oscilloscope records 
the time- variable signal of the dynamic pressure sensor 
(Section l2.5[) . All components are powered via adapters 
fed by a 220 V (50 Hz) supply, secured with a Ground 
Fault Interrupter (GFI), an emergency switch, a master 
fast fuse (5 A), and individual fuses for each component. 



2.2 Pressure- regulated test chamber 

The test chamber (see Figs.[5jD, d) is the heart of the ex- 
periment. It contains demineralized water, where cavi- 



tation bubbles are produced. This test chamber is a cus- 
tomized version of "Vacuum Chamber C" distributed 
by Terra Universal and has two essential properties: 
(1) it is made of acrylic glass, which is optically trans- 
parent as required by the laser-based bubble generator 
and the high-speed camera system; (2) it can hold a 
vacuum, i. e. withstand pressures 100 kPa (1 bar) be- 
low the outside-pressure, as required for a systematic 
study of cavitation bubbles in different pressure envi- 
ronments. The lid of the test chamber can be removed 
to access its inside and exchange the liquid. 

The test chamber is filled to 80% with liquid and to 
20% with air. This 'air' also includes the vapor of the 
liquid at saturation pressure, as well as potential traces 
of laser-generated and bubble-generated gases. The air 
is separated from the liquid by a passive 'gas exchanger' 
(see Fig.[5Ji)- In normal gravity and hyper-gravity (anti- 
parallel to the z-axis), this gas exchanger works as fol- 
lows. The conical shape at the bottom of the exchanger 
forces gas bubbles in the liquid to migrate towards the 
center of the exchanger, from where they escape to the 
air phase across a series of vertical tubes. This perma- 
nent escape channel for gas is particularly important to 
constantly remove traces of laser-generated and bubble- 
generated gases from the liquid. Vice versa, the conical 
shape at the top of the gas exchanger implies that liq- 
uid trapped in the air phase flows down to the center of 
the exchanger, from where it leaks down to the liquid 
through the vertical tubes. The efficiency of this phase 
separation is enhanced by a hydrophobic Teflon coating, 
which prevents small drops and air bubbles to remain 
attached to the gas exchanger by surface tension. In 
microgravity, the small diameter-to-length ratio of the 
vertical tubes (diameter = 4 mm, length = 40 mm) im- 
poses enough friction to prevent the liquid from migrat- 
ing to the air phase under the effect of random gravity- 
fluctuations at the O.Olg-level. This key feature was ver- 
ified experimentally by analyzing the high-speed movies 
and complementary webcam images (Section 12.51) . No 
unwanted liquid/air mixing was found, except for in 
the few cases (< 2%), where the microgravity phase 
was perturbed by flight-turbulences causing an upwards 
gravitational acceleration of a z > O.lg. 

The chamber pressure is lowered below the ambient 
pressure (i. e. the aircraft cabin pressure of 82 kPa) us- 
ing a low-power vacuum pump (Parker BTC-IIS mini- 
pump) . The suction of this pump is connected to the air 
phase of the chamber via a valve. The minimal achiev- 
able air pressure p a ; r in the chamber is 9 kPa. This pres- 
sure is constantly measured by a static pressure sen- 
sor (see Section 12.51) , which covers the whole accessible 
range from 82 kPa down to 9 kPa at an RMS-precision 
of ^0.2 kPa or 1% (whichever is larger). Both the pump 
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and the pressure sensor are connected to the computer, 
allowing the latter to regulate p a i r . Under the effect of 
gravity, the pressure po inside the liquid at the position, 
where the cavitation bubble is generated, differs from 
p a i r due to the weight of the water above the bubble 
center. Explicitly, 



Po = Pair - pa z H, 



(1) 



where p is the density of the liquid, a z is the gravita- 
tional acceleration along the z-axis (a z = — g in normal 
gravity), and H = 70 mm is the height of the water 
between the bubble and the free surface of the liquid. 



2.3 Bubble generation system 

A single cavitation bubble is generated at the center of 
the test chamber by a focused laser pulse. This pulse of 
8 ns duration, 532 nm wavelength, and selectable en- 
ergies up to 230 mJ emerges from a frequency-doubled 
Q-switched Nd:YAG-laser (Quantel CFR 400). We first 
describe the physics behind laser-based bubble genera- 
tors and subsequently elaborate on the particular ad- 
vances incorporated in the present experiment. 

Laser-based techniq ues for generating bubbles have 
been wid ely explored ( LauterbarnL 19721; Lauterborn 
and BolleT'— rr ' ' * ' 



1975 



Tomita and Shimal. Il990l; IVogel et al 



19961: Phi 


iDD and Lauterbornl. Il998l: IVoeel et all. 


1999 


Ohl et all. 


19991: IWolfrum et al 20031: Bruian et al 


2002 



They generally work as follows. A parallel pulsed laser 
beam, typically of a few nano-seconds duration and vis- 
ible color, is focussed inside an optically transparent liq- 
uid medium - say water. Despite the transparency of 
the water, enough energy is absorbed in the focal point 
that a small liquid volume (radius < 0.1 mm) heats 
up to ionizing temperatures (i. e. several 1000 K). This 
liqui d volume hence tr ansform s into a plasma fVoge l 



et al. Il996l : IVoeel et all . Il999l : iBvun and Kwakl . l20o"l 
yielding an enormous pressure of roughly 0.5 GPa per 
1000 K assuming an ideal gas law. The hot plasma ex- 
pands explosively, first by compressing the liquid in 
its immediate environment, thus generating a spheri- 
cal shock wave, then by radially accelerating the liquid. 
The newly formed bubble typically expands to a radius 
i?o 100-times larger than that of the original plasma. 
This enormous expansion (volume factor of 10 6 ), as 
well as the efficient heat transfer through radiation and 
conduction, quickly cool the expanding plasma, which 
hence deionizes and then recondensates. Apart from 
non-condensible gases gene r ated b y the laser pulse (e.g., 
H 2 and O2, see lSato et all . I2013T ). the bubble has 'for- 
gotten' its thermal history by the time it reaches its 
maximal radius. Therefore, a vapor bubble generated 



by a laser pulse is for many practical purposes identical 
to a cavitation bubble, i. e. a bubble that grows due to 
low liquid pressure rather than high gas pressure. 

Since the plasma pressure is per se isotropic, laser- 
generated bubbles expand similarly in all directions, re- 
sulting in spherical bubbles. However, deviations from 
a sphere, can originate from the non-spherical shape of 
the microscopic plasma and from inhomogeneities in the 
liquid pressure. Both sources of asymmetry have been 
minimized in our experiment, as described hereafter. 

Firstly, in order to obtain a highly symmetric laser- 
generated plasma, the parallel laser beam must be fo- 
cussed into a single point at a large angle of conver- 
gence (> 30°), since such an angle reduces the heating 
of water inside the laser beam in the vicinity of the fo- 
cal point. Len ses are not very suitable for such a focus 
(Chapter 4 in iTinguelvl . 120131 ). since they suffer from 
monochromatic aberration at large angles of conver- 
gence and since the quality of the lens-focus depends 
on the wavelength of the laser and on the refractive in- 
dex of the liquid, thus requiring a customized lens for 
each laser-liquid pair. Those issues can be bypassed by 
using a concave parabolic mirror to focus the paral- 
lel laser beam. Therefore, the present experiment uses 
for the first time a parabolic mirror (Edmund Optics, 
2" x 2" 30 deg off-axis parabolic gold mirror, #NT47- 
088), such as shown in Fig. d. This mirror exhibits 
an angle of convergence of about 53°. A gold-coated 
mirror surface was selected to suppress a degradation 
by corrosion. The gold surface absorbs about 30% of the 
laser energy at 532 nm - an acceptable artefact given 
that the surface density of the absorbed energy remains 
a factor 100 below the damage threshold. To avoid that 
a part of the laser beam is reflected back into the laser 
source, which might damage the latter, we chose an off- 
axis mirror as show in Fig.^. The parabolic mirror has 
a declared surface accuracy of A/4 (RMS) and rough- 
ness of less than 17.5 nm (RMS). These high-quality 
specifications allow the laser to be focused in a volume 
smaller than the volume of the generated plasma, hence 
avoiding initial asymmetries of the plasma. 

Secondly, pressure irregularities in the liquid must 
be minimized in the quest for a spherical bubble col- 
lapse. Any pressure gradient Vp in the vicinity of the 
collapsing bubble can cause the formation of a jet di- 
rected against Vp (jObreschkow et a i l2011bh . If the liq- 
uid is at rest before the bubble is generated, then Vp 
can have two origins: a static uniform pressure gradient 
Vp = pa is due to a (gravitational) acceleration a, and 
a dynamic non-uniform Vp is due to the interaction 
between the m oving b ubble surface and nearby bound- 
aries. In fact, iBlakel (1988) showed that bubbles col- 
lapsing near a flat rigid boundary form a jet directed to- 
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wards that boundary, whereas bubbles collapsing near a 
flat free surface form a jet directed away from that sur- 
face. This boundary-induced jet-formation dominates 
over gr avity- induce d jet- formation, if (square of Eq. 8.8 
in Ref. iBlakd . Il988h 



A 



pa 



< 0.2, 



(2) 



RoAp 

where h is the distance between the bubble center and 
the flat boundary (rigid or free surface) , a is the norm of 
the gravitational acceleration, Rq is the maximal bub- 
ble radius, and Ap is the pressure difference between 
Po and the pressure inside the bubble (approximately 
the vapour pressure of the liquid). In most past ex- 
periments A is much smaller than 0.2, either intention- 
ally or due to the difficulty of generating a bubble far 
away from components necessary for the bubble gen- 
eration (e. g. optical lenses). Hence the effects of the 
gravity-induced pressure gradients normally remain ob- 
scured by boundary effects. To decrease these effects, 
such that A > 0.2, the value of h in Eq. ^ must be 
maximized. To do so, we chose a test chamber that 
is large compared to the cavitation bubbles (dimen- 
sions in Fig. (2}i) and a parabolic mirror with a large 
focal distance of 51 mm. Given the off-axis geometry 
of this mirror, the average distance between the focal 
point and the mirror is h = 55 mm. Given water (p — 
10 3 kg m -3 ) in normal gravity (a — g = 9.81 m s~ 2 ), 
Eq. @ then implies that gravity-driven jet-formation 
occurs if RqAp < 150 kg s -2 . In standard units this 
condition for 'gravity-domination' can be rewritten as 
/Ro \ fA£ 
\mm / \bar / 
For example, at a pressure Ap = 0.2 bar all bubbles 
with maximal radii Rq below 7.5 mm will be dominantly 
deformed by gravity, rather than by boundary effects. 

In summary, the generation of a spherical bubble 
requires a parabolic mirror with a large angle of conver- 
gence to maximize the symmetry of the initial plasma 
and large focal length to minimize the effects of nearby 
boundaries. Together these two requirements imply a 
large diameter of the parabolic mirror (here 51 mm). A 
Galilean beam expander (Newport High-Energy Laser 
Beam Expander HB-10X) is used to expand the parallel 
laser beam to this large diameter (see Fig. [5b, c). 

The energy E p of the laser pulse can be adjusted be- 
tween ~ mJ and ~ 200 mJ by varying the time-delay 
between laser-pumping and Q-switching from 500 fxs 
down to 170 /is. Only a minor fraction of this energy 
is transformed into the cavitation bubble, while other 
parts are converted into the shock-wave at the plasma 
formation, converted into heat of the liquid, absorbed 
by the chamber wall and gold mirror, or simply trans- 
mitted across the fluid. The maximal bubble energies 



< 1.5. 



(3) 



E\y that can be obtained with this setup are about 
Et, ~ 12 mJ. These energies are calculated from the 
observed maximal bubble radius via 



E h = ^RlAp, 
where 



Ap=p Q -p v . 



(4) 



(5) 



Po is the water pressure, calculated via Eq. [T] from the 
measured pressure p a j r , and p v is the vapor pressure 
calculated from the measured wa ter temperature via 
Antoine's equation ( Antoinel . 1888 1. 



2.4 High-speed visualization system 

The visualization system contained within the upper 
stage of the rack is schematically represented in Fig.^. 
It consists of a high-speed camera (Photron Fastcam 
SA1.1) fitted with a 135 mm objective (Nikon Zoom- 
Micro 70-180mm f/4.5-5.6D) and exchangeable astro- 
physical filters to analyze the sonoluminescent flash. 
We use a set of three broad-band filters (Baader Anti- 
reflective RGB CCD-filters, 65 x 65 mm) to analyze the 
sonoluminescent flash in three distinct wavebands (R, 
G, B), according to the transmission spectra shown in 
Fig. [3] An additional clear filter is used in order to main- 
tain the same focus when no color filter is applied. This 
RGB-filter system proofed to be useful for an approxi- 
mate characterization of the luminescence temperature. 

The camera operates at inter-frame times down to 
2 /is and exposure times of 370 ns. The spatial resolu- 
tion of the data released with this paper is 69.0 /xm/pixel 
horizontally (y-axis) and 68.7 /xm/pixel vertically (z- 
axis). The typical field-of-view is 17.7 x 17.6 mm 2 (256 x 



l DOW 
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Fig. 3 (Color online) Transmissivity spectra of the optional 
color filters placed in front of the high-speed camera. The 
spectral properties and planarity of these filters satisfy the 
high standards of astronomical photography (Ref.: Baader 
RGBC-CCD 65 x 65 x 3 mm Filters). 
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256 pixels), although it can be expanded to 70.7 x 
70.3 mm 2 (1024 x 1024 pixels). 

A parallel LED background illumination (StageLine 
3W LED-36Spot) provides a white background, against 
which cavitation bubbles and jets appear as black ab- 
sorption features, such as illustrated in Figs.QJi and[2f. 
This background-illumination exhibits the advantage 
that the bubble surface is defined very clearly. Hence, 
the dimensions of the bubbles and jets can be measured 
reliably at a sub-pixel resolutions of 7 /im. Using the 
minimal exposure times of 370 ns, this illumination also 
allows the visualization of shock waves in shadowgraphy 
as shown in Fig. [TJ;. An alternative front illumination 
(not yet available aboard the parabolic flights) , also al- 
lows us to visualize bubbles as bright features in front 
of a dark background. This illumination is more appro- 
priate to study the 3-dimensional morphology and the 
inner structure of the bubbles, such as the thin vertical 
jet visible inside the rebound bubble in Fig. [T]d. 

No illumination is used when studying the sonolu- 
minescent flash from the collapsing bubble itself. The 
CMOS-sensor of the high-speed camera is sufficiently 
sensitive to detect this faint flash on most single bub- 
ble collapses, even when filters are used. Section 13.51 
explains how the signal-to-noise ratio of this data can 
be increased by exploiting the cross-correlation between 
the sonoluminescent pulse and the collapse shock recor- 
ded by the dynamic pressure sensor (Section 12.51) . 

The high-speed camera generates a significant power 
of heat, worth about 100 W, which must be continu- 
ally removed from the closed upper stage of the rack 
through an active air-cooling system. 



2.5 Sensors 

In addition to the high-speed camera an array of sensors 
is used for experimental and security reasons. These 
sensors are described one-by-one in the following. 

2. 5. 1 Sensors for science and control 

• Dynamic pressure sensor: A high-frequency pres- 
sure sensor is used to record the transit of the shock 
waves produced at the initial plasma generation and 
at the subsequent collapses of the principal cavita- 
tion bubble and the cascade of rebound bubbles. 
Our first experiments with already published re- 
sults (jObreschkow et all . l2011bl . I2OI2I : iTineuelv et al 
2012[ ) use a custom-designed piezo-resi stive pressure 



sensor described in Section 4.2.1 of (Hasm atuchi 



to accurately measure the timing of the shock pas- 
sage, but the shock front cannot be sampled in de- 
tail. Measurements of the shock energy therefore re- 
quire extrapolation techniqu es, which proofed to b e 
useful but rather uncertain ( Tinguelv et al 2012 '). 
To improve these measurements the pressure sen- 
sor was replaced by a modern hydrophone (1.0 mm 
needle with 28 /xm sensor, manufactured by Pre- 
cision Acoustics PAL). The bandwidth of this hy- 
drophone extends to 12 MHz, hence allowing a de- 
tailed sampling of the shock front, although the ac- 
tual peak-pressure, typically on the order of hun- 
dreds of MPa, remains inaccessible. The pressure 
signal, represented by an electrical voltage is pre- 
amplified and recorded by an oscilloscope (LeCroy 
WaveRunner 6050A, 500 MHz bandwidth). 
Static pressure sensor: A gas pressure sensor (Phid- 
gets Gas Pressure Sensor, #1115), connected to the 
chamber-side of the vacuum pump (Fig. [5fc) , moni- 
tors the pressure of the air inside the test chamber at 
an RMS-precision of 0.2 kPa or 1% of the measured 
pressure (whichever is larger). Using the non-linear 
calibration displayed in Fig. |4j the nominal linear 
pressure range of (20-250) kPa, could be extended 
down to 9 kPa without loss of precision. For a mea- 
sured static pressure pphidgets the true air pressure 
is calculated by 

PPhidgcts 



Pair = 



' PPhidgcts \ 

v 11.5 kPaJ 



(6) 



201 if ). This sensor has a baseband bandwidth (con- 



stant gain to 25 kHz, resonant at 100 kHz) and a 
pressure range (linear in 0-500 kPa), which suffice 



Chamber temperature sensor: A thermometer is used 
to measure the water temperature, from which the 
vapor pressure is calculated. Our first experim ents 
( Obreschkow et a]| . l2011b1 ; lTinguelv et~all2012h used 
a sensor (Phidgcts Precision Temperature Sensor, 
#1124) attached to the outside wall of the test cham- 
ber (Fig. This sensor has a statistical precision 
of 0.5 K, but since it is not in direct contact with 
the water, the measurements are subject to system- 
atic uncertainties and time-delays. To improve the 
measurements the sensor was subsequently replaced 
by a water-resistant thermistor (USB Thermistor 
DTU6024C-004-S), inserted into the water near a 
corner of the test chamber. This sensor reliably mea- 
sures the water temperature at a statistical accuracy 
of 0.1 K. 

Accelerometer: A 3-axis accelcromctcr (Phidgets Ac- 
celerometer 3- Axis, #1059), fixed on the rack, re- 
cords the inertial acceleration a (including gravity) 
in the range [-3g, +3g] = [-29.4 ms~ 2 , +29.4 ms~ 2 ] 
in each dimension, at an RMS-precision of 0.0019<? 
(a x and a z ) and 0.0029(7 (a y ). The bandwith of 30 Hz 
lies safely above the fastest g-jitter of the aircraft 
around 10 Hz. This accelerometer is used both to 
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Measurement of static pressure sensor Pphidgets W 

Fig. 4 (Color online) Non-linear calibration used to extend 
the nominal pressure range of the static pressure sensor (20- 
250 kPa) down to 9 kPa. The data points show the true air 
pressure p a i r applied to the static pressure sensor against the 
measurements pphidgets of the latter. The non-linear function 
(solid line) of Eq. ([6]) provides an excellent fit to these data. 



trigger the generation of a cavitation bubble at the 
desired level of gravity and for the post-processing, 
if we require a gravity-measurement rather than a 
nominal value associated with the particular flight 
maneuver. An example of a measurement of dz as a 
function of time is shown in Section 12.61 (Fig. [5]) . 
Webcam: A USB-webcam (Logitech Quickcam Pro, 
960 x 720 pixels, 30 frames per second) is installed 
inside the upper stage of the experimental rack in 
order to safely monitor the setup, while the exper- 
iment is running. All major components of the up- 
per stage are visible on the webcam image, such as 
shown in Fig. [SJ 

Manual trigger: A hand-hold button, stowed in the 
lower stage of the rack, can be used to release an ex- 
perimental cycle (bubble generation and high-speed 
imaging). This button serves as a backup in case of 
an in-flight failure of the automatic trigger. 



2.5.2 Sensors for experimental security 




Fig. 5 (Color online) Webcam image allowing a safe moni- 
toring of the experiment when the box is closed during laser 
operation. The components described in Fig.[5]that are visible 
in the Webcam image are labelled. 



directly to the computer, prevents the latter from 
commanding the laser, if the lid is open. 
Light sensors: Two light sensors (Phidgets Precision 
Light Sensor #1127) monitor the illuminance inside 
the upper rack stage at the 1 lux-level. If traces 
of light are still measured, while the magnetic lid 
sensor is closed and the LED illumination system is 
turned off, the laser is stopped automatically. This 
mechanism provides an additional level of security, 
excluding the possibility that the laser is fired if the 
upper rack stage is damaged in such away that light 
can escape even though the lid is closed. 
Laser temperature sensors: A temperature sensor 
(Phidgets Precision Temperature Sensor, #1124) at- 
tached to the laser head (Fig. [2};) and four ther- 
mometers built into the laser-box by the manufac- 
turer control the laser temperature at all times. If 
any of these sensors indicates a dangerous temper- 
ature, the laser is stopped immediately. 
Amperemeter: An inductive amperemeter (Phidgets 
30 Amp Current Sensor AC/DC, #1122) monitors 
the current drawn by the experiment at 1.5% ac- 
curacy. Any current irregularities below the 5 A 
threshold of the master fuse, are displayed as warn- 
ings on the computer. 



• Mechanical lid sensors: Two electro-mechanical con- 
tact sensors, directly connected to the "interlock" 
connector of the laser power unit, ensure that no 
laser pulse can ever be fired while the lid of the up- 
per rack stage is open. 

• Magnetic lid sensor: As an additional security mea- 
sure, a magnetic contact sensor (Phidgets Minia- 
ture Magnetic Contact Switch BR-2023), connected 



2.6 Parabolic flights 

The experiment is performed on parabolic flights flown 
by the aircraft "A300 zero-g" - so far during the 53rd 
(year 2010) and 56th (year 2012) parabolic flight cam- 
paign of the European Space Agency (ESA). The typi- 
cal layout of a flight campaign is as follows. Over three 
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consecutive flight days, the aircraft performs 3 x 3f = 93 
parabolas, offering about 20 s of microgravity each. De- 
pending on demand, the steep turns are added to gen- 
erate stable levels of hyper-gravity at 1.2g, lAg, 1.6g, 
and I.85. A single parabola and a steep turn at 1-Sg 
are illustrated in Fig. [5] together with a representative 
measurement of the vertical acceleration, recorded by 
the accelerometer. 

The parabolas are flown in a ballistic manner, i. e. at 
a constant horizontal speed, thus imitating the center- 
of-gravity of motion of a free-falling body in the ab- 
sence of air drag. By doing so, the aircraft experiences 
weightlessness in the same way as a spacecraft in orbit. 
There are, however, small gravity-perturbations due to 
tiny deviations in the aircraft's trajectory. This residual 
g-jitter exhibits an RMS-amplitude of O.OI5 at frequen- 
cies of 1-10 Hz. Furthermore, the slow pitching of the 
aircraft causes a small centrifugal acceleration that is 
strongest at the nose and at the tail of the aircraft, 
where it reaches about O.Olg. Our experiment has an 
allocated space near the aircraft's center of gravity to 
avoid this source of systematic perturbation. 



2.7 Trigger sequence 

Each experimental cycle with a single cavitation bubble 
follows the same sequence. A software trigger activates 
the laser unit as soon as the user-defined level of gravity 
is reached and stabilized for 1 s. This unit initiates the 
pumping of the laser and uses an internal precision clock 
to actuate the Q-switch at the preset time-delay. When 
the Q-switch is activated, thus releasing the laser pulse, 
a signal simultaneously triggers the high-speed camera 
and the oscilloscope recording the dynamic pressure. 
We chose to synchronize the camera and oscilloscope 
with the Q-switch to ensure that they are always syn- 
chronized with the cavitation bubble, independently of 
the Q-switch delay. 



3.1 Parameter space overview 

A total of 4887 single cavitation bubbles were generated 
during the flight maneuvers, using demineralized water 
at room temperature. Of those bubbles, 247 were ob- 
served with background illumination, while the remain- 
ing 4640 bubbles were recorded without illumination to 
detect the sonoluminescent flash (details in Section l3~5j) . 
Three parameters could be varied independently: 

• The energy of the laser was varied to modulate the 
bubble energy Ey, between 1.1 mJ and 11.9 mJ. 

• The air pressure p a j r of the chamber was adjusted 
to vary Ap between 6.5 kPa and 81.2 kPa. 

• The gravitational acceleration a z slowly alters be- 
tween Og and 1.8g during the nights, causing the 
pressure gradient | Vp| to vary between and 1.8p<7. 

As shown in Section 13.41 the evolution of our bub- 
blcs is closely described by Rayleigh's theory (jRavleighl . 
1917J). Therefore, given E\, and Ap, the maximal bub- 
ble radius Rq (1.5 — 7.0 mm) and its collapse time T c 
(157 — 2439 /is) are uniquely determined via Eq. (Q| and 



To = fR 




(7) 



where the Rayleigh-factor is / « 0.91468. Using Eqs. (j4j) 
and ||ZJ), any two of the variables {E^, Ap, R ,T C } fully 
determine the remaining two. In logarithmic space these 
relations are linear, i.e., 



log-Bb = ki 
logT c = k 2 



3 log Ro + log Ap 
logi? - 0.51ogZ\p 



(8) 
(9) 



with constants k\ = log(47r/3) and k 2 = log/+0.51ogp. 
Using these relations we can represent all four parame- 
ters {£"bj Ap, Rq, T c } in a two-dimensional logarithmic 
plot as shown in Fig. [71 The situations probed by our 
experiment are represented as points, whose color indi- 
cates the level of gravity. 



3.2 Overview of a bubble's life 



3 Results 

This section describes the data gathered specifically 
during the 53rd parabolic flight campaign conducted 
in October 2010. We first overview the sampled param- 
eter space ( Section l3.ip and illustrate the life-cycle of 
a single cavitation bubble (Section 13.21) . Thereafter, we 
successively address technical details associated with 
different stages in this life-cycle (Sections l3.3H3.5l) . Note 
that sy stematic scientific analyses are published sepa- 
rately (jObreschkow et all . l2011bl . l2012UTinguelv et all . 
20121 . and forthcoming). 



Fig. [5] illustrates the time- resolved measurements char- 
acterizing the evolution of a single bubble. 

Fig. shows the high-speed movie of a bubble spec- 
ified in the caption of Fig. [51 This bubble is generated 
at the time t — /is. At t = 2 /is the "primary" shock, 
driven by the expanding plasma, can be seen about 
3 mm from the bubble center. This shock is visible 
due to the density-depende nce of the water's refrac- 
tive index (|Cho et at 120011 ). The apparent thickness 
of the shock of 0.5 ± 0.1 mm can be explained by the 
smearing of a sound wave (c ~ 1480 ms _1 ) traveling for 
the exposure time of 370 ns. Hence, the intrinsic shock 
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Fig. 7 (Color online) Overview of the bubble parameters 
probed on parabolic flights. The data is binned in three grav- 
ity bins (colors) according to the nominal gravity level of the 
flight manoeuvre. The four parameters E^, Ap, Rq, and T c 
are related via Eqs. @ and (0) as demonstrated in Section 
13.41 Therefore those four parameters can be reduced to only 
two dimensions, as done in this plot. 
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Fig. 9 Interleaved high-speed sequence of the bubble gen- 
eration from a hot, laser-generated plasma. The laser pulse 
at t = lasts for 8 ns. (Bubble energy E^ w 3 mj, driving 
pressure Ap ?s 20 kPa, exposure times are 370 ns) 



sure signal also indicates significant shocks at the first 
and second bubble collapse. In addition to these clear 
features in the pressure signal, the signal also shows 
various spurious oscillations, partially associated with 
higher-order reflections, partially caused by the sensor's 
eigenmodes around and beyond 100 kHz. 

The different stages of the bubble evolution will now 
be discussed in detail over the Sections 13.31 to [ 



3.3 Phase I: Bubble generation 



thickness must be smaller than 0.1 mm. On the interval 
/is — 2741 /zs the bubble grows and collapses. At the 
collapse point, a "secondary" shock wave is emitted. 
The collapse is followed by a rebound, during which a 
gravity-driven jet emerges from the bubble. 

Fig. \E]p shows a synthetic 3-color movie based on 
2400 individual movies of bubbles similar to that shown 
in Fig.|5Ji, but filmed without background-illumination. 
Those movies were taken while placing the three color 
filters (R, G, B) successively in front of the camera, such 
that each color band yields about 800 movies. Those 
movies, all taken at inter-frame spacings and exposure 
times of 66.7 /is (15 kHz), were then rescaled in time 
to match the bubble collapse time of Fig. [Sh. and av- 
eraged to obtain a noise-reduced signal. The resulting 
color movie clearly shows the white sonoluminescent 
radiation at the collapse point, which indicates a tem- 
perature on the order of 10 4 K in the black body ap- 
proximation. The beginning of the movie shows a bright 
green spot, which is an inevitable after-glow of the in- 
tense green (532 nm) laser pulse. 

Fig. |8j; shows the dynamic pressure signal. At f w 
47 /Lis, the plasma-shock appears on the dynamic pres- 
sure signal, the delay being due to the distance of 70 mm 
between the shock center and the sensor (Fig. [5Ji) . At 
t w 142 /is, the pressure exhibits a negative peak, asso- 
ciated with the stratification wave generated when the 
shock is reflected at the free water surface. The pres- 



The bubbles are generated by a 8 ns laser pulse focussed 
inside the liquid (Section l2.3p . Fig. [5] shows three inter- 
leaved photographs of the initial bubble growth at an 
inter-frame spacing of only 270 ns. These photographs 
were extracted from three different movies of three sim- 
ilar bubbles (bubble energy w 3 mj, driving pres- 
sure Ap w 20 kPa) . The exact timing of the images was 
determined from the timing of the primary shock on 
the dynamic pressure sensor. In this way it is possible 
to mimic inter-frame spacings much shorter than the 
minimal inter-frame spacing of 1.5 /is. 

Fig. unveils that a spherical shock wave is gen- 
erated while the laser-heated region shines at visible 
wavelengths. The brightness of this light-emitting re- 
gion causes the camera sensor to saturate and blur. The 
white disk at the center of Fig. [5^ is therefore much 
larger than the actual light-emitting region. The exis- 
tence of radiation at visible wavelengths witnesses tem- 
peratures of several 1000 K, hence the light-emitting 
medium must be ionized. On Fig. \§]p this plasma has 
cooled down to temperatures that still emit visible light, 
without, however, saturating the camera sensor. At this 
stage, a dark proto-bubble centered about the plasma 
becomes visible. Stated differently, Fig.[9j) demonstrates 
that there is a moment in the bubble-generation pro- 
cess, where the bubble has already formed and the shock 
has already detached from the bubble, although the 
bubble gas is still hot enough to be partially ionized. 
However, only 270 ns later, i. e. at about 0.1% of the 
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Fig. 10 (Color online) Growth and collapse motion of four 
distinct cavitation bubbles, represented in normalized coordi- 
nates. The solid line is the prediction of the Rayleigh model. 
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Fig. 11 (Color online) Comparison of the half-life times Tj/2 
of 163 cavitation bubbles measured using the dynamic pres- 
sure sensor (T-P^J 5 ) and the high-speed movie (T™° v ) , respec- 
tively. The close agreement between these independent mea- 
surements demonstrates the small measurement uncertainties 
of both types of measurements. 



Rayleigh time (Eq. [7]), the plasma has fully recombined 
(Fig. |9j;) , leaving a bubble fille d with wa t er vap or and 
minor amounts of other gases (jSato et all . I2013T ). 



normalized to the maximal radius Ro and the total life- 
time To. Although these data correspond to four dif- 
ferent experimental conditions, the normalized evolu- 
tions are congruent and closely approximated by the 
Rayleigh model (solid line in Fig. 1101) . In fact, we have 
recently shown that the collapse motion of some of our 
bubbles fits the R ayleigh model within les s than 0.1% 
in terms of radius (| Obreschkow et all . |2012| ) . This close 
agreement with the Rayleigh model suggests that the 
effects of surface tension and viscosity are negligible - 
a feature attributed to the large (Ro > 1 mm) size of 
the bubbles considered in this work. 

Fig. [TU] also reveals that R(t) is nearly symmetric in 
time. Hence the collapse time T c , i.e, the duration from 
R = Ro to R ~ 0, is well approximated by the half- 
life time T1/2 = To/2. Slight asymmetries in the time- 
evolution are physically interesting, but lie beyond the 
technical description envisaged in this paper. Instead, 
we shall now show that the experiment can measure 
T1/2 with a relative uncertainty below 10 -3 . In general, 
T X j2 can be measured independently from the pressure 
signals and high-speed movies. The pressure-based ap- 
proach measures the time-difference between the pri- 
mary shock, emitted at the bubble generation, and the 
secondary shock, emitted at the final stage of the bub- 
ble collapse. Since those two shocks imprint a similar 
response on the sensor (see Fig. [Bt), their separation in 
time can be very accurately determined using an auto- 
correlation applied to the spline-interpolated voltage- 
signal. Alternatively, I\/2 can also be measured directly 
from the high-speed movies. These movies use different 
inter-frame spacings from 1.5 p& to 20 fj,s. By measuring 
the bubble radii on each frame at sub-pixel resolution 
and applying a Rayleigh-interpolation to estimate the 
continuous time-evolution of the radius, Ti/ 2 can be 
measured with an uncertainty much smaller than the 
inter- frame times. 

We can then compare the half-life times T^° s and 
T™° v extracted from the pressure signals and movies, 
respectively. To this end we use all 163 cavitation bub- 
bles that have high-speed movies and full pressure data. 
These bubbles cover half- life times of T1/2 = 157 — 
2439 /is, energies of -Eh = 1 — 10 mJ, and pressures of 
Ap = 9 — 80 kPa. The comparison (Fig. lTTj) reveals that 
^l/T anc ^ ^i/2 V systematically agree within a standard 
deviation of 



1 l/2 1 l/2 



380 ns. 



(10) 



3.4 Phase II: Bubble growth and collapse 

Once generated, the bubbles grow and collapse. Fig. ITUl 
shows this evolution in terms of the bubble radius R(t), 



Assuming that the measurement uncertainties of T^J° S 
and T™° v are uncorrelated, this result implies that Ty^ 
and T™° v both have a statistical measurement uncer- 
tainty below 380 ns. In other words, both the pressure 
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Fig. 12 (Color online) Measured Rayleigh factors /, defined 
in Eq. [7] compared against the theoretical value. 



sensor and the high-speed camera measure n at a rel- 
ative precision between 0.01% and 0.1%, depending on 
the actual value of Ti/2- This is an extreme accuracy 
compared to the nominal time-resolution of the pres- 
sure sensor (~ 10 /is at resonance) and the inter-frame 
spacing (1.5 — 20 fis) of the movies, thus demonstrat- 
ing the power of using smart time-interpolations. Since 
those interpolations are simpler in the case of T^, 2 S than 
T™° v , we hereafter set Ti/ 2 — T^° s , whenever a reliable 
pressure signal is available, and Ti/ 2 — Tw™ otherwise. 

Upon assuming T c — T\/2 and provided parallel 
measurements of Rq from the high-speed movies, we 
can then check if the pairs {T c , Rq} are systematically 
consistent with the Rayleigh model. To do so, we solve 
Eq. for the Rayleigh factor / while propagating the 
measurement uncertainties of T c , Rq, and Ap. The val- 
ues of / hence obtained are consistent with the theoret- 
ical value (Fig.ll2|). Gravity has no obvious effect on the 
life-time of the cavitation bubbles besides the indirect 
effect of changing Ap via Eq. Jl} . 



3.5 Phase III: Inflection point and rebound 

The collapse of a cavitation bubble progresses until a vi- 
olent increase in the gas pressure abruptly decelerates 
the collapse motion (jAkhatov et all . l200lh . The com- 
pres sion of the non-con densabl e bubble gas f Tinguelv 
et al, 20121 ) increases the internal energy of this gas to 
ionizing temperatures, thus causing a sonoluminescent 
flash. Simultaneously, the immense gas pressure signif- 
icantly compresses the surrounding liquid, thus driving 
the secondary shock away from the bubble. Finally, the 
gas pressure stops the collapse motion and inverts it, 
thus causing a rebound bubble. As outlined in the intro- 
duction, the key objective of our setup is to analyze the 
energy partition between spherical rebound, micro-jet, 
shock wave, and thermal effects such as luminescence 
(see Fig.[T]). To achieve this goal, we record the evolving 
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Fig. 13 Collapse and rebound of two cavitation bubbles 
imaged against a bright background. The top row shows a 
bubble in microgravity, where the rebound remains spheri- 
cal, while the bottom row shows a bubble in normal gravity, 
where the rebound is deformed b y a jet propagating again st 
the gravity vector, [adopted from Obreschkow et al (2011b)] 



shape of the rebound, the pressure of the shock wave, 
and the sonoluminescent flash. 

The evolving rebound bubble is imaged by the grey- 
scale high-speed camera used with background illumi- 
nation. Those visualizations reveal a gravity-driven jet 
(Fig. I X3[) , which disappears in microgravity conditions 
and gets enhanced in hyper-gravity. This jet is in actual 
fact a vapor envelop of a thin liquid micro-jet, visible 
in Fig. [lb. These jets are caused by the gravity- induced 
pressure gradient Vp = pa. However, since the jets can- 
not 'know' that the pressure gradient is due to gravity, 
we can consider gravity as a generic tool to study the 
effects of any uniform pressure gradient Vp. Further- 
more, Vp is the first factor in the spatial Taylor expan- 
sion of any inhomogeneous pressure field. Hence, this 
study also applies, to first order, to any inhomogeneous 
pressure field. We recently presented detailed investiga- 
tions of the jet a s a function of Vp and othe r experimen- 
tal parameters ( Obreschkow et aj 2011b ). This study 
found that the volume of the vapor jet, normalized to 
the volume of the spherical component of the rebound, 
depends linearly on the non-dimensional parameter 



C = \V P \R /Ap. 



(11) 



A theoretical model confirmed this scaling law. Ground- 
based follow-up experiments with viscous mixtures of 
water and glycerine further showed that the jet-size is 
independent of the viscosity of the liquid. The jet's inde- 
pendence of viscosity, as well as its independence of sur- 
face tension can be seen as a consequence of the global 
conservation law of momentum (Kelvin-impulse) dur- 
ing the inflection of the collapse motion f Obreschkow 
et al,|2011bj). A" different explanation for the insignifi- 
cance of viscosity and surface tension is that their ef- 
fects are masked by inertial forces, which increase more 
rapidly as R(t) — » 0, according to the Rayleigh-Plesset 
equation (jPlesset and Chapman! . 1 19711 ). 
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The shock wave is registered on the dynamic pres- 
sure sensor with a time-delay of 47 [is. Although the 
piezo-resistive sensor used on the 53rd parabolic flight 
campaign is unable to sample the wavefront, the inte- 
grated response of this sensor, J dt s 2 (i), can be used 
to estimate the energ y of the shock. In recent work 
(jTinguelv et all . |2012| ) , we calibrated this method and 
analyzed systematic variations of the energy ratio be- 
tween the shock wave and the rebound bubble. This 
empirical study, backed-up with a first-order model for 
liquid shocks, suggests that the energy partition be- 
tween shock and rebound is fully determined by the 
single non-dimensional parameter 



7 6 Z\p 



1/7 
Pg 



(pc 2 ) 1 



-1/7 



(12) 



where p g is the pressure of the non-condensable gas 
at the maximal bubble radius, 7 ps 1.3 is the adiabatic 
index of the non-condensable gas, p is the liquid density, 
and c is the liquid's speed of sound. An example of two 
bubbles with different values of £ is shown in Fig. [T4l 

Finally, the greyscale high-speed camera can cap- 
ture the sonoluminescent flash instead of imaging the 
rebound bubble, if the illumination is turned off. The 
luminescence only ever appears on a single frame of 
the movie. Hence this luminescence lasts much shorter 
than the minimal inter-frame spacing of about 1.5 fj,s, in 
agreement with detailed theoretical and experimental 
studies suggesting typ i cal du rations of less than 10 ns 
( Lauterborn and Kurz . 2010h . Our setup can image the 
luminescence either directly, through any of the three 
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Fig. 14 (Color online) Two bubbles with different normal- 
ized rebound radii. The fraction of the initial bubble en- 
ergy recovered in the rebound bubble is related to the non- 
dimensional parameter £, defined in Eq. (112ft . These two bub- 
bles have values of £ K, 600 (black *) and £ ~ 70 (red +), the 
difference originating from different values of Ap. The other 
bubble parameters are given in the caption of Fig. 1101 



color filters (RGB, see Fig. or through a new spec- 
tral analyzer not yet available on the 53rd parabolic 
flight campaign. Only one filter can be used at a time; 
hence three-color studies must rely on several bubbles 
produced under identical conditions and observed suc- 
cessively in the R, G, and B band. The problem with 
this method is the limited reproducibility of the lumi- 
nescence. In fact, we found that bubbles of identical en- 
ergy and driving pressure Ap, thus identical radius and 
life-time, vary in their luminescence brightness about 
10-times more than expected from pure shot noise con- 
siderations. These brightness fluctuations are probably 
related to the microscopic size of the sonoluminescent 
plasma, which makes it highly sensitive to minor per- 
turbations and easily obscured by nuclei and impuri- 
ties in the water. Any robust result must therefore rely 
on statistical averages of many single bubbles. For this 
reason, 20-times more experiments were conducted to 
measure luminescence than to image the rebound bub- 
ble on the 53rd parabolic flight campaign. Correspond- 
ing scientific results will be presented in forthcoming 
work. 



4 Online data 

Most data presented in this paper can be accessed on- 
line at bubbles. epH.ch/data. So far, this online data 
contains 247 cavitation bubbles imaged against a white 
background during the 53rd parabolic flight campaign. 
These bubbles have been numbered from 'cavityOOOOl' 
to l cavity00247' in chronological order. Each bubble has 
a single zip-file, containing up to eight data files: a log- 
file with header data, a file with the dynamic pressure 
data, and six files for the high-speed movies. Some bub- 
bles are missing either the pressure file of the movie files 
due to experimental difficulties in-flight. 



4.1 Log-files (cavityOOOOOJog.txt) 

The log-files are ascii-files describing the experimental 
conditions of the bubble in all available detail. An ex- 
ample of a log-file is given in Tab. [TJ In particular, the 
log-files contain all available measurements of the static 
sensors. Furthermore, they contain comments on prob- 
lems with certain sensors, as well as potential comments 
made by the experimenter. At the end of each log-file, 
we list basic post-processing results, including the max- 
imal bubble radius Rq, the collapse time T c , as well as 
the driving pressure Ap. Given those values and their 
measurement uncertainties, we also provide the most 
likely values of Rq , T c , and Ap assuming that they sat- 
isfy the Rayleigh model in Eq. (JTJ). Those values, forced 
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to satisfy Eq. ([7]), are called reference values. Explicitly, 
they have been obtained by maximizing the total prob- 
ability function p(Ro)p(T c = JRq \J p / Ap) p(Ap) , where 
p(Rq), p(T c ), and p(Ap) denote the individual proba- 
bility distributions associated with the individual mea- 
surements of Ro, T c , and Ap, respectively. Most exper- 
imental values are provided with uncertainty intervals. 
They represent standard deviations where the noise is 
considered Gaussian, and smallest intervals containing 
67% of the events, where the noise is non-Gaussian. 



4.2 Pressure files (cavity00000_pressure.txt) 

The pressure files are ascii-files containing the pre-amp- 
lified voltages of the piezo-resistive dynamic pressor 
sensor, sampled at 2.5 MHz using an oscilloscope (Le- 
Croy WaveRunner 6050A, 500 MHz bandwidth). Each 
file contains 20,000 data points, hence covering a time- 
interval of 8 ms, sufficient to capture the primary and 
secondary shock of all bubbles. The data is organized 
in two columns. The first column gives the time in s 
relative to the instant of the laser pulse that generates 
the bubble. The shock impacts the pressure sensor with 
a time-delay of about 47 p,s, due to the 70 mm distance 
between the bubble center and the pressure sensor. The 
second column gives the electrical voltages recorded by 
the oscilloscope. These voltages have been calibrated to 
13.55 /iV Pa -1 . However, the sensor's response is too 
slow and too limited in dynamic range to measure the 
actual shock pressure. The pressure signal can nonethe- 
less be used to measure the life-time of the bubble and 
to estimate the shock energy as outlined in Section [3J 
Some bubbles exhibit dynamic pressure data with- 
out clear shocks and/or clear low- frequency oscillations 
on the order of 1 kHz. Those data are affected by spu- 
rious air bubbles sitting on top of the pressure sensi- 
tive piezo-element. Typically, the data deteriorated in 
this way yield values below 5 in the 'Dynamic pressure 
signal-to-noise' entry of the log-files. 



4.3 Movie files (cavity00000_movieX) 

There are up to six movie files per cavity specified 
by the placeholder 'X'. They are three actual video- 
files and three corresponding header-files, readable by 
the software Photron Fastcam Viewer (PFV), provided 
by the manufacturer of the high-speed camera. The 
three files are the original 12-bit greyscale video (cav- 
ity00000_moviel2bit.mraw), readable by PFV, as well 
as two standard 8-bit greyscale movies in avi-format, 
readable by most video software. The first avi-file (cav- 
ity00000_movie8bit.avi) is a simple resampling of the 



original 12-bit scale. For most applications it is vir- 
tually identical to the original 12-bit movie, but gen- 
erally more easy to process. The second avi-file (cav- 
ity00000_movie8bit_enhanced.avi) is graphically enhan- 
ced movie with a flattened background brightness and 
reduced noise achieved via dark frame subtraction. 



5 Summary and future prospects 

In this work, we have presented an experimental setup 
able to produce highly spherical cavitation bubbles in 
conditions so rigorously isotropic that the pressure gra- 
dient of normal gravity is the most important source of 
asymmetry. In fact, the experiment can clearly visual- 
ize a gravity-driven jet, propagating against the vector 
of gravitational acceleration. In the past, such observa- 
tions of gravity-driven jets could only be obtained in the 
case of large bubbles (Rn > 1 cm) collapsing in hyper- 
gravity ([Benjamin and EllisL l 19661 ) , hence requiring less 
spherical initial conditions. Given that normal gravity 
is the most important anisotropy of our bubbles, we 
then conducted the whole experiment in microgravity 
aboard parabolic flights. This setting thus allowed us 
to generate bubbles more spherical than theoretically 
possible in any ground-based experiment. 

As argued in the introduction, this experiment of- 
fers an ideal platform to uncover a general model for 
the energy partition between rebound, jet, shock, and 
thermal effects during the collapse of a single cavitation 
bubble. Such a model would make a crucial tool for 
reducing cavitation damage and optimizing beneficial 
applications of cavitation. While this paper is focussed 
mostly on technical aspects, we have also illustrated two 
important scientific results: The strength of the micro- 
jet is propor tional to the non-dimensi onal parameter £ 
of Eq. JIT} (|Obreschkow et all l2011bl ). and the energy 
partition between the rebound and the shock depends 
monotoni cally on the non-dim ensional parameter £ of 
Eq. (H2) (jTineuelv et all . l2012h . 

In making some data obtained during the 53rd para- 
bolic flight campaign (October 2010) available online, 
we open a possibility for other research groups to work 
with exclusive microgravity data. Moreover, next gen- 
eration data obtained on the 56th parabolic flight cam- 
paign (May 2012) and forthcoming ones will allow a 
more detailed view of the luminescence. Ultimately, this 
research is aimed at unifying the quickly diversifying 
field of cavitation research in providing a global model 
for the energy distribution between all cavitation re- 
lated effects. 
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Fig. 2 (Color online) Overview of the experimental setup installed inside the Airbus A300 zero-g (53rd ESA parabolic flight 
campaign), (a) Photo of the entire mechanical structure ("rack") fixed inside the aircraft. The lower stage of the rack hosts 
electronic control/acquisition systems and power supplies. The upper stage of the rack consists of a closed box containing the 
test chamber, the laser-system for cavity generation, cameras and other sensors, (b) Photo of the subsystems located inside 
the upper stage of the rack, (c) Schematic top-view of the upper stage of the rack, (d) Side-view of the vacuum vessel, which 
constitutes the test chamber for the cavitation bubbles, (e) Schema of the aircraft showing the location of the experiment, (f) 
Image of the cavitation bubble as seen by the high-speed camera. (ALL) All six panels show the three basis vectors (x,y,z) of 
the Cartesian coordinate system of the aircraft. By definition, this basis is 'left-handed'. Note, however, that on the camera 
image (e), the coordinates appear 'right-handed' due to the mirror between the test chamber and the camera objective. 
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Steep turn 




Fig. 6 (Color online) Illustration of a parabola and a steep turn flown by the aircraft A300 zero-g. A representative measure- 
ment of the vertical acceleration — a z as a function of time is plotted underneath the flight trajectory. 
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Fig. 8 (Color online) Illustration of the time-resolved measurements characterizing the evolution of single cavitation bubbles. 
Fig. [8^ shows selected 370 ns-exposures of a high-speed movie of a bubble evolving in a static pressure Ap = 8.98 ± 0.03 kPa 
and a gravitational acceleration \a\ = (1.847 ± 0.005)g. This bubble attains a maximal radius Ro = 4.549 ± 0.007 mm and 
exhibits a collapse time of T c = 1410.83_ 'q2 fJS. The background-illumination used in this visualization turns the bubble and 
shocks into dark absorption features. Fig.[8p shows a color- movie obtained without illumination, in order to visualize the light 
emitted by the bubble itself (luminescence). Each color channel (R, G, B) of this movie is the average of about 800 single 
movies obtained by using the color filters (Fig. on 2400 single cavitation bubbles. Fig. [BJ; displays the voltage recorded by 
the dynamic piezo-pressure sensor (Section 12. 5|l during the evolution of the bubble shown in Fig. [5Ji. The physical phenomena 
revealed in these figures are explained in Section 13.21 
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Name 



Additional explanation 



Value 



Unique cavity ID 
Version of data reduction 

Date of cavity generation 

Time of cavity generation 
[UTC+2h] 

Flight maneuver 



Flight maneuver sub-index 

Problems found in data inspec- 
tion 

High-speed movie available 
Control sensors available 

Dynamic pressure available 
Dynamic pressure signal-to-noise 

Pulse rate 



Pulse index 



Dynamic pressure mode 



Q-switch delay [us] 



Nominal vessel pressure [kPa] 



Unique 5-digit identifier of the bubble; the bubbles presented 
here range from 00001 to 04887. 

Version of the code used to reduce the movie- and pressure- 
files and to produce the log-files. 



'dayX' shows the flight day X = 1,2,3; 'paraYY' means 
that the bubble was produced in modified gravity during the 
parabola YY = 00, .., 30 or in the normal gravity succeeding 
that parabola; 'turnZZ' means that the bubble was produced 
in the hyper-gravity of the steep turn ZZ = 01, ..,08. 

Bubbles belonging to the same 'Flight maneuver' are distin- 
guished with an increasing sub-index, ranging from 1 to 200. 



Is 'false' if the movie is not available due to technical problems 
in-flight. 

Is 'false' if data from the control sensors (e. g. acceleration, 
static pressure, temperatures) is not available due to technical 
problems in-flight. 

Is 'false' if data from the dynamic pressure sensor is not avail- 
able due to technical problems in-flight. 

Ratio between the peak voltage and the RMS-noise of the 
dynamic pressure sensor. Values below 10 indicate that the 
pressure signal is deteriorated by gas on the piezo-sensor. 

'single' means that two successive bubbles are spaced by sev- 
eral seconds, such that any interference between those bub- 
bles can be ignored; '7Hz' means that bubbles were generated 
at a 7 Hz rate, hence gases produced by a bubble may not be 
completely removed until the next bubble is formed. 

Equals '1' if the 'Pulse rate' is 'single'; and counts from '1' 
upwards if the 'Pulse rate' is '7Hz'. 

If 'DC, the pressure signal is recorded without filter, if 'AC, 
the pressure signal has been processed by a high-pass filter. 
Note that the offset tension is negative in 'DC and zero in 
'AC. 

Delay between laser pumping and Q-switching. Larger values 
indicate less laser energy. Range: 170 (corresponding to ~ 
200 mj pulses) to 500 (vanishing pulse energy). 

Target value of the pressure p a ir, which the computer tries 
to achieve; only use this value, if no value is available for 
'Measured vessel pressure'. 



Measured vessel pressure [kPa] Pressure p a i r measured by the static pressure sensor. 



Derived pressure at cavity level 
[kPa] 

Vapor pressure of the water [kPa] 



Pressure po at the bubble center before bubble generation, 
obtained via Eq. fljl. 

Vapor pressure p v derived from a measurement of the water 
temperature at the beginning of each flight day. 



Measured ambient pressure [kPa] Pressure in the aircraft cabin. 

Nominal gravity |g| [9.8m/s"2] Norm of the gravitational acceleration that the aircraft tried 

to achieve during the respective flight maneuver. Only use 
this values if no 'Measured gravity' is available. 

Measured gravity |g| [9.8m/s~2] Measured norm a = \a\ of the gravitational acceleration. 



00147 
2.0 

21-Oct-2010 
09:59:13.995 

day3turn04 



true 



true 



true 



18.2 



single 



275 
10.3 

10.32±0.15 
11.58±0.15 

2.46 

83.00±0.5 
1.8 

1.847±0.005 
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Name Additional explanation Value 



Measured gravity g_x [9.8m/s"2] 


Measured gravitational acceleration a x . 


-0.194+0.005 


Measured gravity g_y [9.8m/s"2] 


Measured gravitational acceleration a y . 


-0.017±0.005 


Measured gravity g_z [9.8m/s~2] 


Measured gravitational acceleration a z . 


-1.837±0.005 


Movie size 


Size in pixels of the frames of the high-speed movie. 


iODXZOD 


Number of frames 


Number of frames in the high-speed movie. 


600 


Frame rate [Hz] 


Frame rate of the high-speed movie. 


67500 


Shutter speed [s] 


Exposure time of each frame. 


1/2700000 


Microns per pixel (horizontal) 


Horizontal size of each pixel in the high-speed movie. 


69.9 


Microns per pixel (vertical) 


Vertical size of each pixel in the high-speed movie. 


69.6 


LED backlight 


Is 'on', if the back-illumination is used (e. g. Fig. |8ji); is 'off' 
of no illumination is used (e. g. Fig. [Hp). 


Oil 


Camera diaphragm 


F-number of the camera objective, that is the ratio between 
the aperture diameter and the focal length. 


11 


Camera filter 


Type of filter placed in front of the camera as shown in Fig. [2};; 
'clear', if a clear filter is used that corresponds to no filter, 
but conserves the focus; 'red', 'green', or 'blue' if one of the 
three color filters described in Fig. [3] is used. 


clear 


rrame 01 luminescence 


If 'LED backlight' is 'off', this is the frame number in the 
high-speed movie, on which the sonoluminescent light is ei- 
ther clearly detected or expected from the collapse shock mea- 
sured by the dynamic pressure sensor. 


IN /A 


Sonoluminescent signal [summed 
pixel value] 


Number of counts of the sonoluminescent signal. WARNING: 
It is necessary to check in the high-speed movies whether in- 
dividual pixels are saturated (i. e. at value 4192 or 1.0 de- 
pending on the movie- type). 


IN /A 


Measured cavity radius [mm] 


Maximal radius _Ro of the cavitation bubble as measured from 
the high-speed movie. 


4.548±0.007 


Measured half-life time [us] 


Best measurement of the half-life time T1/2 (half the time 
from the bubble generation to the first collapse). This time 
is measured either from the dynamic pressure sensor or the 
high-speed movie, if no dynamic pressure is available. 


1370.40±0.10 


Measured collapse time [us] 


Estimated collapse time T c from the maximal bubble radius 
to the collapse point. This time is obtained from T1/2 via 
T c /s = 1.107(Ti /2/s) 1 ' , an equation calibrated on a subset 
of 20 cavitation bubbles imaged with increase spatial and 
temporal resolution. 


1410.91±0.10 


Source of measured collapse time 


(see above) 


Dynamic pressure sen- 
sor 


Measured delta pressure [kPa] 


Pressure Ap calculated via Eq. lf5l) . 


ci 1 o_l_n 1 ^ 
y.lz±U.10 


Reference cavity radius [mm] 


Most likely value of the maximal bubble radius _Ro , given the 
measurements of _Ro, T c , and Ap, and imposing Eq. (Q). 


4.551+0.007-0.007 


Reference collapse time [us] 


Most likely value of the collapse time T c , given the measure- 
ments of Rq, T c , and Ap, and imposing Eq. Q. 


1410.81+0.20-0.00 


Reference delta pressure [kPa] 


Most likely value of the driving pressure Ap, given the mea- 
surements of _Ro, T c , and Ap, and imposing Eq. ([?)). 


8.70+0.03-0.03 


Reference energy of cavitation 
bubble [mj] 


Value i?b calculated via Eq. Q , given the reference values of 
_Ro and Ap. 


3.44+0.03-0.03 



Table 1: Explanation of the log-files; example of cavity00147dog.txt. 
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